VEmtTMm  Of 

mnasfTY 

STAWFORO.  eAtJFGRNIA 


RELATIONSHIPS  AMONG  CLASSES  OF 
SPHERICAL  MATRIX  DISTRIBUTIONS 


TECHNICAL  REPORT  NO.  10 


KAI-TAI  FANG  AND  HAN-FENG  CHEN 
INSTITUTE  OF  APPLIED  MATHEMATICS,  ACADEMIA  SINICA 

AND 

WUHAN  UNIVERSITY 


APRIL  1984 


U.S.  ARMY  RESEARCH  OFFICE 
CONTRACT  DAAG  29-82-K-0156 


DEPARTMENT  OF  STATISTICS 
STANFORD  UNIVERSITY 
STANFORD,  CALIFORNIA 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


RELATIONSHIPS  AMONG  CLASSES  OF  SPHERICAL  MATRIX  DISTRIBUTIONS 


TECHNICAL  REPORT  NO.  10 


KAI-TAI  FANG  AND  HAN-FENG  CHEN 
INSTITUTE  OF  APPLIED  MATHEMATICS,  ACADEMIA  SINICA 

AND 

WUHAN  UNIVERSITY 


APRIL  1984 


U.S.  ARMY  RESEARCH  OFFICE 
CONTRACT  DAAG  29-82-K-0156 


DEPARTMENT  OF  STATISTICS 
STANFORD  UNIVERSITY 
STANFORD,  CALIFORNIA 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


THE  VIEv'J,  OPINIONS,  AND/OR  FINDINGS  CONTAINED  IN  THIS  REPORT  ARE  THOSE  OF 
THE  AUTHOR(S)  AND  SHOULD  NOT  BE  CONSTRUED  AS  AN  OFFICIAL  DEPARTMENT  OF  THE 
ARMY  POSITION,  POLICY,  OR  DECISION,  UNLESS  SO  DESIGNATED  BY  OTHER  DOCUMEN¬ 
TATION. 


1. 


Introduction. 


The  theory  of  multivariate  analysis  has  been  based  mainly  on  the 
normal  population.  Statisticians  have  been  trying  to  extend  the  sample 
theory  in  multivariate  analysis  to  cases  of  the  observations  being 
not  necessarily  either  normal  or  independent.  In  the  last  decade, 
es:pecially  the  last  five  years,  many  statisticians  have  been 
interested  in  a  specific  class  of  distributions,  one  of  elliptically 
contoured  distributions,  and  found  that  this  class  has  many  properties 
similar  to  the  normal  distribution. 

If  the  c.f.  (charcteristic  function)  of  an  n-dimensional  random 
vector  X  has  the  form  exp(it’]i)(()(t'Zt) ,  where  y:  nxl,  E:  nxn 
and  Z  ^  0,  we  say  that  x  has  an  elliptically  contoured  distribution 
with  parameters  y,  Z,  and  ip,  and  write  x  ^  EC  (y,Z,(})).  When  y  =  0 

5  ~  5n’  call  EC^(0,I^,(j))  a  spherical  distribution  and  write 
X  S  ((})),  because  x  -v  S  ((}))  iff  x  =  Fx  for  each  T  e  0(n),  where 
0(n)  is  the  set  of  nxn  orthogonal  matrices  and  the  notation  "x  =  y" 
means  that  x  and  y  have  the  same  distribution. 

As  an  extension  of  multivariate  normal  sampling  theory,  several 
classes  of  spherical  matrix  distributions  are  defined  and  have  been 
discussed  by  many  authors.  Here  are  the  main  three  classes. 

Definition  1.  Let  X  =  (x^, .  . .  ,x^)  be  an  nxp  random  matrix. 

(1)  ^2.  ^  ^  5  every  T  e  0(n)}.  The  X  in  is  called 

left-spherical  by  Dawid  (1977). 

(2)  ^2~  ^5’  ^El^l’ ' * ■ ’^p^p^  ^  Cx^,...,x^)  for  every  e  0(n), 

1  =  1, . . . ,p}. 

(3)  =  {X:  r(vecX)=  vec  X  for  every  V  e  0(np)}.  Each  of  F^, 
end  F^  is  called  a  class  of  spherical  matrix  distributions. 
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These  classes  contain  many  important  distributions,  such  as  the 
multivariate  normal  distributions,  the  multivariate  t-distributions , 


the  multivariate  Beta-distributions  and  the  multivariate  stable  lav7s. 

The  class  was  studied  by  Dawid  (1977,  1978),  Fraser  and  Ng  (1980), 

Jensen  and  Good  (1981),  and  Kariya  (1981a  and  b).  The  class  F2  was 
defined  by  Anderson  and  Fang  (1982b) ,  and  the  class  F^  was  discussed 
by  Chmielewski  (1980),  Kariya  (1981a),  Jensen  and  Good  (1981)  and 
Anderson  and  Fang  (1982b  and  c).  They  have  found  that  many  statistics 
are  invariant  in  these  classes.  However  the  relationships  among  these 
classes  are  not  yet  very  clear.  Therefore  it  may  be  valuable  to 
consider  the  relationships  among  them. 

Throughout  this  paper,  X  =  (x^,...,Xp)  =  (x^^^ , . . . ,x^^^ ) '  =  (x^^) 

denotes  an  nxp  random  matrix  with  n  ^  p,  denotes  the  nxn 

identity  matrix,  diag  (a^, . . .  ,a^)  denotes  an  n^n  diagonal  matrix  with 

diagonal  elements  a,,..., a  ,  A’,  rkA  and  trA  denote  the  transpose  of 

i  n  ~  ~  ~ 

A,  the  rank  of  A  and  the  trace  of  A,  u^’^^  denotes  a  random  vector 
which  is  uniformly  distributed  on  the  unit  sphere  in  R^,  and  (t't) 
denotes  its  c.f. 

In  Section  2  some  basic  properties  about  them  are  listed.  Sections 
3  and  4  are  the  main  part  of  the  paper.  In  the  last  section  we  stimmarize 
the  invariant  distributions  in  these  classes  as  a  table  for  applications. 


2.  Preliminary. 

In  this  section  we  recall  some  basic  properties  of  F^,  F^,  and  F^ 
which  will-be  used  frequently  in  this  paper.  The  following  lemmas  and 
Theorem  1  are  from  Dawid  (1977),  Kariya  (1981),  and  Anderson  and  Fang  (1982b). 
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Let 


matrix. 


T  -  (t^,,..,t  )  be  an  nxp 
~  ~P 

Lemma  1. 


(1) 

X  £  F- 

"v  X 

iff 

the 

c .  f 

.  of 

X 

has  the  form 

(|)(T’T); 

(2) 

5  ^  ^2 

iff 

the 

c .  f 

.  of 

X 

has  the  form 

ct>(tjti,.. 

(3) 

X  £  F^ 

iff 

the 

c.f 

.  of 

X 

has  the  form 

(j)(trT'T) . 

Lemma  2 

.  Suppose 

X 

has  a  pdf  (probability  density 

function) . 

Then 

(1) 

X  £  F^ 

iff 

the 

pdf 

of 

X 

has  the  form 

f(X’X); 

(2) 

X  £  F^ 

iff 

the 

pdf 

of 

X 

has  the  form 

f  (xj^x^, ... 

,x'x  ) ; 

(3) 

X  £  F 

iff 

the 

pdf 

of 

X 

has  the  form 

f (trX’X) . 

(a) 


Theorem  1. 

J  ^  if  X  =  U^A,  where  A  and  are  independent, 


and  U'U,  =  I  ; 

~1~1  ~p 

(b)  X  e  F  iff 

~  Zt 

are  independent, 


(c)  X  e  F^  iff 


vec 


d 


d 

5  "  H2~’  5  =  diag(R^, , . . ,R^)  and  U2 =  (u^, 

>  0>  i  =  l,...,p,  and  u  , ...,u  i.i.d,  = 
1  —  ~p 

X  =  XU^,  where  R  ^  0  and  are  independent, 


The  distribution  of  is  called  the  uniform  distribution.  The 

matrices  U^,  U2  and  play  roles  of  coordinate  systems  in  F^, 

F2  and  F^,  respectively.  In  this  paper  when  we  write  X  =  XU2, 

5  ~  H2~’  5  ~  Hl~’  have  the  meaning  in  Theorem  1  unless  we 

make  another  explanation. 
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Lemma  3.  The  distribution  of  X  is  fully  determined  by  that  of 


X'X  if  X  £  Fj,  by  that  of  (xix  , )  if  X  £  F2,  and  by  that 
of  trX'X  if  X  £  F^. 

Lemma  3  shows  us  that  if  X  £  F^,  Y  £  F^,  and  X'X  =  Y'Y,  then 
X  =  Y.  Similar  statements  hold  for  ^2  ^2' 

The  following  properties  of  the  operation  "="  are  given  by  Anderson 
and  Fang  (1982a) : 

(1)  If  X  =  Y  and  f^(*)»  i  =  l,...,m,  are  Borel  functions,  then 
(f^(X),...,f^(X))  =  (f^(p,...,f^(Y)). 

(2)  If  z  is  independent  of  X  and  Y,  respectively,  then 

(a)  X  =  Y  implies  zX  =  zY; 

(b)  if  P(z  >  0)  =  1  and  the  c.f.  of  log  z  ^ 

d  d 

almost  all  t  then  zX  =  zY  implies  X  =  Y. 


Lemma 


_4.  Suppose  X  £  F^  and  P(lx'xl=0)  =  0,  then 


(1)  X  =  U^A  =  U^B  implies  A  =  B,  where  A  and  B  are  two  upper 
triangular  matrices  with  positive  diagonal  elements; 

(2)  X  =  QT,  where  Q'Q  =  1^  and  T  is  an  upper  triangular  matrix  with 
positive  diagonal  elements,  implies  Q  =  U-  and  Q  is  independent  of  T. 


Proof.  The  existence  of  A  in  (1)  follows  from  the  argument  of 
Dawid  (1977).  We  consider  mapping  f:  A  f (A)  =  A'A,  where  A  is  an 
upper  triangular  matrix  with  positive  diagonal  elements.  Clearly,  f 
is  a  one-to-one  mapping.  As  a'’A  =  B’B,  then  (1)  follows  from 
E[h(A)]  =  E[h(f"^(A'A))]  =  E[h(f“^(B'B))]  =  E[h(B)]  for  each  Borel 
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function  h  _>  0.  Note  that  if  |x'x|  ^  0,  there  is  a  unique  decom¬ 
position  X  =  QT,  where  Q'Q  =  and  T  is  an  upper  triangular 
matrix  with  positive  diagonal  elements.  Let  the  function  g(X)  =  (Q,T), 
we  have  (Q,p  =  g(X)  =  g(Uj^A)  =  (U^,A)  as  X  =  U^A,  which  completes 
the  proof.  Q.E.D. 

From  now  on  when  we  write  X  =  U^A  for  X  £  F^,  we  always  consider 
that  A  is  an  upper  triangular  matrix  with  nonnegative  diagonal  elements. 

3 .  Relationships  among  ,  F^ ,  F^ . 

By  Definition  1,  clearly  =  ^2  ^  ^3* 
differences  among  them  at  all?  In  this  section,  we  will  discuss  them 
in  the  following  aspects:  the  coordinate  system  and  the  coordinate 
transformations,  marginal  distributions,  marginal  densities  and  spheri¬ 
city.  First  of  all,  start  with  U^,  U2  and  U^. 

Lemma  5.  Suppose  X  e  F^  and  P(x^=0)  =  0,  i=l,...,p.  Then 
X  e  f 2  X  satisfies  the  following  conditions: 

(1)  x^/|lxj|  ,...,Xp/||Xp|l  are  independent;  and 

^Il5lll  ’'"'’lljpll  ^  >  •  •  •  >^p/ ll^pll  )  3.re  Independent. 

Proof.  The  assertion  follows  from  Anderson  and  Fang  (1982b).  Q.E.D. 

Corollary  1.  i  F2. 

Proof .  As  =  I^  and  e  F^,  =  (u^,...,u^)  where 

Uij*..,Up  are  not  independent,  the  corollary  follows  from  Lemma  5.  Q.E.D. 
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Lemma  6.  U2  i- 

Proof .  Suppose  U2  =  , .  . .  )  '  eF^,  then 

U2  =  RU^  for  some  R  ^  0  being  independent  of  U^.  As  ,u^ 

are  independent  and  u..,  (i  =  l,...,n)  has  a  spherical  distribution, 

~  viy 

the  distribution  of  must  be  normal  and  the  distribution  of  U2 

must  be  normal  by  Kelker  (1970).  The  contradiction  proves  the 
theorem.  Q.E.D. 


The  following  example  shows  us  that  the  condition  (2)  in  Lemma  5 
is  necessary,  and  its  distribution  belongs  to  Fj^,  but  not  to  T- 


Example  1.  Let  X  be  an  nxp  random  matrix  with  a  pdf  f (X)  = 
where  c  is  a  constant.  We  want  to  prove  that 

p/ll;pli )  =  H 

(R^,...,R^)  and 

U  are  independent,  then  f (X)  =  f (UR)  can  be  written  in  the  form  of 
f(X)  =  g^(R)g2(U);  similarly  ll^+X'xl  =  h^(R)h2(U)  for  some  two 
-*  Ul.i.VA  Ha*  Take  R 

1  2  ~  ~p  ^  ^  ^p- 


c|l  +X'xl“^^^'^P^ 

,  where  c 

~P  ~  ' 

(llxjl  ,.  .  .,  llxpll  ) 

E  (R^,...,Rp) 

not  independent. 

It  is  easy  to 

functions  h^  and  h„.  Take  R  =  I  and  U'U  =  I  ,  respectively. 


to  show 


I  +X'xl  =  ll  +RU’Ur1  =  k|l  +R  Ml  +u'u|  , 


for  each  R  =  diag(R^, 
and  ulu,  =  1,  i  =  1,. 
t  >  0,  we  have 


. .  ,R  )  with  R-  >0,...,R  >0, 
P  1  ’  ’  p  ’ 

.,p,  where  k  is  a  constant. 


U  —  (u^  , .  . .  jU,  ) , 
~1  ~P 

Let  R  =  t  , 
~P 
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k 1 1  I  +I  I  1 1  +U 'U I  =  I t^I  +U’U I ,  for  each  t  >  0, 

~p  ~p  ~p  ~  ~  ~p  ~  ^ 

i.  e. 

k(l+t^)P  n  (1+A.)  =  n  (t^+X.),  for  each  t  >  0, 
i=l  ^  i=l  ^ 

where  are  the  eigenvalues  of  U'U.  Obviously,  this  is 

impossible.  This  contradiction  reveals  that  R  and  U  are  not  inde¬ 
pendent  . 

3.1  Coordinate  transformations. 

In  this  paragraph  we  try  to  give  the  stochastic  representations  for 
?!»  ^2  ^3  under  the  same  "coordinate  system".  Maybe  they  can  help 

us  to  understand  these  classes  more  clearly.  The  generalized  Dirichlet 
distribution  defined  by  Anderson  and  Fang  (1982a)  will  be  used.  Consider 

J  O 

a  random  vector  (z^  , .  . .  ,z  )  *  such  that  (z,  , .  .  .  ,z  )  S  R'^(ci  ...  d  )  , 

where  O^R'VFCO,  R  is  independent  of  (d^ ,  •  .  .  >  d^) ,  d^  +  •  •  •  +  d^  =  1 , 

and  (d-,...,d  ^ D  (aT  , . .  .  ,a  ;a  )  (Dirichlet  distribution)  with 

1  m-1  m  1  m-1  m 

a^,...,a^  >  0  and  n  =  2(a^+*--+a^)  is  an  integer;  we  write 

(zi, .  .  .  ,Zm_i)  *'^l’-"’^m^  G^(a^, .  . .  ,a^;F)  . 

Anderson  and  Fang  (1982b)  pointed  out  that  if  =  U2R  with 

R  =  diag(R  ,  .  .  .  ,R  )  ,  then  (R^, .  . .  ,R^  ,  )  '^^  D  (-^, .  . .  ,-^;^)  and 
~  ±  p  i  p-1  p  2  2  2 

R  +...+R^  =  l.  Also  X  =  RU  e  F„  and  R%F(*)  iff  X  S  u.R 

with  R^  =  diag(R^, - ,R^)  'v  g  (-^j  . . .  ,-|n;F)  . 

Let  Y  =  (y^j)  =  (y-L»'**jyp)  be  an  nxp  random  matrix  with 
i.i.d.  elements  and  y^^  'u  N(0,1).  There  exists  an  upper  triangular 
matrix  T  =  (t..)  with  positive  diagonal  elements  such  that  T'T  =  Y'Y. 
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This  is  the  famous  Bartlett  decomposition.  It  is  a  well-known  fact 


that  independent,  t^. ^  N(0,1)  for  i  <  j 

and  Xn_i+i’  1  =  Let  t^  =  (t^^  , . .  .  ,t^^  )  ’  , 

j  =  1 , . . . , p .  Then  we  have 


Theorem  2.  Suppose  =  U^A,  where  A  is  an  upper  triangular 
matrix  with  positive  diagonal  elements.  Then  A  =  TR  where 
R  =  diagdi  t^ll  , . .  .  ,  lltpll  ). 

Proof.  As  (y]_/l|yilU •  •  •  »yp/||ypll)  "  H2  ^  Hi~  I'l  " 


have  yly  =  t't  ,  1  =  l,...,p,  and 


u^A  =  YR  ^  ^  ^ 


where  QT  is  the  Bartlett  decomposition  for  (y^^ . y^)  and  Q'Q  - 

By  Lemma  4,  we  have  TR~^  =  A  which  completes  the  proof.  Q.E.D. 


Remark.  Let  A  =  (a..),  a.  =  (a, . , . . . ,a. . ) ' ,  a.  =  (a. a.  ,  .)' 
-  ~  11  ~i  li  ’  11  ’  ~i  li  1-1,1 

(2)  2  2 

and  a^  =  (a^^, . . . ,a^^) ,  i  =  l,...,p.  Then  by  Theorem  2,  we  obtain  the 
following  facts: 

(1)  a^,...,a  are  independent; 

1  ^p 


,(n) 


*  d 

^k  =  i^k’ 

where 

and 

ik 


(3)  a^^^  Dj^Cy,  ...,-|-,^(n-k)),  k  =  2,...,p. 

Corollary  1.  X  e  ^2  iff  X  =  U^AR,  where  U^,  A,  and 
R  =  diag(R^, . . . ,Rp)  ^  0  are  independent,  and  A  is  given  by  Theorem  2. 

Theorem  3.  Let  U3  =  U^B;  then  B  =  T/(trT'p^^^. 
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i  H 


Proof.  By  using  the  above  notation,  we  have  Y*(trY*Y)^^^  = 
and  T'T/trT'T  =  Y'Y/trY'Y  =  B’B.  Note  B  is  also  an  upper  triangular 
matrix  with  positive  diagonal  elements,  then  B  =  T/ (trT Q.E.D. 

Corollary  1.  X  e  F„  iff  X  =  RU^B,  where  R  >  0,  U, ,  and  B 
are  independent,  and  B  is  given  by  Theorem  3. 

3 . 2  Classes  of  marginal  distributions. 

Let  (i=  1,2,3)  denote  the  set  of  first  columns  of  X's  in 

F  (i= 1,2,3),  i.e.,  x  £  F^  iff  there  exist  x„,...,x  such  that 

5  ~  ’^2’ '  * ' ’5p^  ^  ^i  (i~L>2,3).  Similarly,  F^  indicates  a  set  of 

the  first  row  vector  of  X  in  F.  (i=  1,2,3).  Clearly  F^  =  F^  o  F^ 

and  =  F^  3  F^,  since  F^  ^  ^2  ^  ^3*  ^2  "  ^1’ 

X  £  F^  let  X2,...,Xp  be  p-1  nxl  random  vectors  such  that 

J’52’'*'’2p  iid;  thus  X  =  (x,X2, . . . ,x^)  e  F2.  Here  we  use  a 

useful  fact  that  if  X  =  (x  ,...,x  )  with  Xt,...,x  i.i.d.  and 

~i  -'P  ~1  ~P 

x^  't,  EC^(0,  I^,(})) ,  then  X  £  ¥2'  Naturally,  one  may  ask  if  F^  =  F^ 
holds.  But  that  is  not  true.  First,  we  need  the  following  lemma. 

Lemma  7 .  (t't)  with  t  £  R^  ^  is  not  an  n+l-dimensional  c.f. 

n  ~  ~  -V/ 

Proof .  Suppose  fi^(t’t)  with  t  £  R^"^^  is  an  (n+l)-dimensional 
c.f.  Then  there  exists  a  distribution  function  F  such  that  (c.f. 
Cambanis,  Huang  and  Simons  (1981)) 

(3.1)  fi^(u)  = 


l^_^^(ur  )  dF(r),  u  ^ 
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i.e. , 


=  Ru  ,  where  R  is  independent  of  u  ,  R  >  0,  and  u 
~n  ~n  —  ’  ~n 

is  the  subvector  of  with  the  first  n  components.  Since 

has  a  pdf  and  P(R=0)  =  P(Ru  =0)  =  P(u^"^=0)  =  0,  has  a  pdf. 

»S»X1  ~ 

This  is  a  contradiction,  which  completes  the  proof.  Q.E.D. 

c  c 

Theorem  4.  The  set  is  a  proper  subset  of  F2  if  p  >  1. 

Proof .  Let  u  =  u^’^^  Clearly,  u  £  F^,  we  want  to  point  out 

u  i.  F^.  Suppose  u  £  F^;  then  there  exist  U2,...,u^  such  that 

U  =  (u,u„,...,u  )  £  F„.  Let  ({)(trT'T)  denote  the  c.f.  of  U.  Then 

-v/p  J  ^  ^  ~ 

the  c.f.  of  u  is  (})(t’t-)  =  9,  (t't,),  t.  £  r’^,  i.e.,  (])(*)  =  (•). 

~  n  ~JL~1  n 

That  means  that  (t't)  =  ({)(t't),  t  £  R^^,  is  a  c.f.  By  Lemma  7, 

n  -V  -s/  ^ 

the  contradiction  proves  the  theorem.  Q.E.D. 


Let  us  consider  the  rox<?  marginal  distributions.  First,  we  want  to 

T'X*  T'X’ 

point  out  that  is  a  proper  subset  of  F^  . 

Lemma  8.  Suppose  X  =  (x^,...,x^)  =  (x^^  , .  .  .  ,x^^^ )  '  £  F2  and  the 
covariance  of  exists.  Then 

(1)  Cov(x^^^  ,x^^  ^ where  is  a  diagonal  matrix  and 

6 .  .  =  1 ,  and  6..  =0,  i^j,  i,j=l,...,n,  and 
1 1 

2  2 

(2)  Cov(x.,x.)  =  6.. a.. I  ,  where  a.,  will  be  given  in  the  proof, 

~i’~j  ij  ii~n’  11  t-  > 

i,3  =  l,...,p. 


Proof .  Clearly,  x^^,...,x^^^  are  identically  distributed  and 


(3.2) 


X  =  U2R  =  (R^u^, . . . ,RpU^)  , 
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where  R,  are  independent.  As  EU2  =  0,  we  have  =  0, 


and  Ex  =0  for  k =  1, . . . ,n;  j  =  1, . . , ,p.  By  (3.2) 


The  first  assertion  follows  from  Eu  =  0  for  i  j ;  k  =  l,...,p. 

±K  J  K 
2 

Similarly,  Ex  x'  =  ER.R.Eu.u!  =  6.,ER.I  /n.  and  the  Lemma  follows. 

1  2  ij  i~n 

Eet  x^j^^,.,.,x^^^  be  i.i.d.,  x^^<^  ^  N(0,E)  and  S  is  no t  a 


(1) 


diagonal  matrix;  by  Lemma  8,  then  x^^^  i  F^,  but  x^^^  e  F^.  Thus 


:(l) 


F2  is  a  proper  subset  of  F^. 


Theorem  5.  Suppose  X  =  (x^^  , . .  .  ,x^^^ )  ’  e  then  X  C  F^  iff 


5(1)  ®  ^3- 


Proof.  The  "only  if"  part  is  obvious.  Suppose  £  ^3?  then 


(1) 

•+5(n)£(n)>  “  * 


x^j  has  a  c.f.  where  <{>(t^^t^y 

c.f.  in  r'^^.  On  the  other  hand,  since  X  £  f 2,  X  has  a  c.f. 
i^(t' t  , .  . .  t' t  ),  Hence  we  must  have 

~1~1  ~p~p 


(3.3)  (j)(r^+- •  .+rp)  =  tp  (r^, .  . .  ,rp)  ,  for  r^  ^  0,  i  =  l,...,p. 


because  they  are  all  the  c.f.  of  (3-3),  we  have 

cl)(t]^t^+---+t^tp)  =  ii;(tj^t^,.  ..,t^tp)  for  all  t^  £  r’^,  i  =  l,. 
X  £  F^.  The  theorem  follows.  Q.E.D. 


,p,  i.e. 


Corollary  1.  Suppose  X  £  F2.  Then  £  P3  iff  Sp((j)) 


Proof .  The  "only  if"  part  is  trivial.  Now  suppose  S^(<}>)' 

By  Theorem  5,  we  get  X  £  F^  since  X  £  F2  and  x^^  £  F^.  Q.E.D. 
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Corollary  2.  The  first  row  is  not  in  r^. 

Proof .  Assume  ^  ^2’  then  there  exists  Y  such  that 

X  =  (u^^,Y)'  e  A.S  is  right  spherical  (Dawid  (1977)),  there- 

fore  ^(1)  ^  ^3  ^  ^  ^2  Corollary  1  of 

Theorem  5.  However,  it  is  impossible  (cf.  the  following  example). 
Hence,  ^  ^2*  Q-^-D. 


Example  2.  Suppose  X  =  U^A,  where  and  A  are  independent, 

nd  A  =  diag(a^, .  .  .  ,ap)  ,  0  <  p^  =  P(a^=l)  =  1-P(a^=0)  <  1,  i  =  l,...,p, 
^+* • •+ap  =  1.  Clearly  X  £  F^,  but  X  ^  F^  by  Theorem  2.  We,  however. 


want  to  show  “  '^^(1)  ^  ^2’  where  is.  the  first  row  of 

2  2 

It  can  be  shown  that  the  c.f.  of  u.-,.  is  Q,  (ti  +  ’ • ‘H-t  )  by  the 

~(1)  n^  1  P  ^ 


sphericity  of  And  ^ 


(3.4) 


2x 


Q  (a?t?+. .•+a^t^)  dF(a, ,...,a  ) 
nil  p  p  1  p 


1=1 


P.fi  (tp 

1  XI  1 


By  (3.4),  we  have 


,n 


ij;(t't  , .  . .  ,t' t  )  =  I  p  (tit  ),  t.  e  R  ,  i  =  l,...,p 

A.  L/  LI  ,  ^  i  Lx  ~  ^  ~  JL  ^ 

^  ^  1=1 


As  (tit.)  is  a  c.f.  in  r'^  and  E?  ,  p .  =  1,  p .  >  0,  i  =  l,...,p, 
n  ~i~i  1=1  *^^1  *^1 


hence  i{i(tl  t., , . .  . ,  t' t  )  is  the  c.f.  of  some  Y  in  F„  and  y 


(1)  ~(i)’ 


where  Y^)  is  the  first  row  of  Y;  that  means  £  F^. 


(1) 
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Example  2  shows  us  that  F^,  related  to  F^,  cannot  be  charac¬ 
terized  by  its  row  marginal  distributions.  But  for  related  to 

F^,  it  can  by  Theorem  5.  Further,  it  is  easy  to  show  that  if 

5’  ?  ^  ^2  ^(1)  ~  ~(1)’  have  X  =  Y.  However,  there  is  no 

such  property  for  F^. 


3 . 3  Marginal  densities. 

Let  XeF^,  1=1,2,3.  In  general  it  is  not  necessary  that  X 

has  a  density.  If  X  satisfies  some  suitable  additional  condition  on 

X,  it  will  have  marginal  densities. 

Suppose  X  =  RU^  £  F^;  if  P(X=0)  =  P(R=0)  =  0,  then  all 

marginal  densities  exist  (Kelker  (1970)).  Suppose  X  =  (R,u,  ,...,R  u  )£F^: 

l--i  P'^P  2 

if  P(x.=0)  =  P(R^=0)  =0,  x^  has  all  marginal  densities;  if 

~  Oj  then  X  has  marginal  densities  of  a  set  of 

elements  such  that  at  least  one  element  in  each  colum  of  X  has  been 
deleted.  Also,  we  can  prove  that  if  X  £  F  and  P(|x'x|=0)  =  0,  then 

(^ll’'--’\-l,l’Xi2’"'’^n-2,2”'*’^lp’*'*’^n-p,p^  and  all  its  subsets 
have  marginal  densities. 

Further,  if  X  =  U2R  z  F 2,  it  is  easy  to  see  that  X  has  a  pdf 

^X ^5l5l ’  *  '  ■  ’ 5pjp ^  ^  there  exists  the 

following  relationship  between  them  (cf.  Zhang  and  Fang  (1982),  Ch.  9): 


Similarly,  suppose  X  -  U^A  £  F^,  where  A  =  (^-^^ )  is  an  upper  trian¬ 
gular  matrix  with  positive  diagonal  elements;  then  X  has  a  pdf  f  (X'X) 
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iff  A  has  a  pdf  f^(A),  and  is  related  to 


£^(A)  .  2V''/2-p(p-1)/4  I  ^ 

i=i  5  -  i=i 


(cf.  Srivastava  and  Khatri  (1979)). 


as  follows: 
r  (•|■(n-i+l)  ). 


4.  The  Class  of  Spherical  Matrix  Distributions. 

Let  A>^0  be  a  pxp  matrix  and  X-  >  •  •  ^  >  X  >0  be  the 

1  —  —  p  — 

eigenvalues  of  A.  We  write  X(A)  =  diag(X^,. . . .X^) . 

Definition  2.  (Dawid  (1977)).  A  random  nxp  matrix  X  is  called 
spherical  if  X  and  X'  are  left-spherical,  i.e.,  PXQ  =  X  for  each 
P£0(n)  and  Q  e  0(p).  Denote  =  {X:  X  is  spherical}. 


Lemma  9.  (Dawid).  X  £  F  iff  X  =  U  AV,  where  U, ,A,  and  V 

[ 

“P 


are  independent,  A  =  X( (X'X)^'^^) ,  TV  =  V,  V'V  =  I  for  each  F  E  0(p). 


The  class  of  F  was  studied  by  Dawid.  The  c.f.  of  X  in  F 
a  ~  s 

must  have  the  form  ({)(X(T'T)),  because  the  maximum  invariant  of 

T(T:nXp)  under  the  transformation  PTQ  for  each  P  £  0(n)  and  each. 

Q  £  0(p),  is  X(T'T).  It  is  easy  to  see  F,  o  F  =  F_  and  U,  E  F  , 
~  ~~  ls3~ls 

If  X  £  F_,  it  is  not  necessary  that  X  £  F  ,  and  vice  versa.  For 
^  s 

example,  the  X  in  Example  1  belongs  to  Fg,  but  X  i  ^2’ 

following  theorem,  we  see  that  U„  ^  F  . 

s 


Theorem  6.  =  f  n  F„. 

-  J  s  2 

Proof .  Clearly,  ^3  ^2'  Conversely,  if  X  £  F^  n  F2,  the 
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fact  X  -  (^(1)  >  •  •  •  >  '  £  Fg  implies  n.  S^((f))  and  X  e 

from  Corollary  1  of  Theorem  5.  The  theorem  follows.  Q.E.D. 


Theorem  7 .  Let 
meaning  in  Lemma  9. 
then 


=  U^AV,  where  U^,A,  and  V  have  the 
2 

Denote  A  =  diag  (A^  , . . ,  ,  A.  >  X  >0, 

~  X  p  X  —  —  p  — 


,Wp)/(w^+---+Wp), 


where  Wj^j...,Wp  are  p  eigenvalues  of  W  W^CnjI^),  and 
(A^,...,A  has  a  joint  density 


(3.5)  -J^hsssm - . 

n  r(^)r(^) 

a=0  ^  ^ 


Kn-p-1) 


n  (A. -A.) 


lli <  j iP 


1  2 


(A.  >  •••  >  A  ,  >  0  and  A  =  1-A  -•••-A  ,  >  0)  . 

1  p-1  p  1  p-1 


and  (A  , ...,A  ,)  is  independent  of  w  =  w,+'‘’4w  . 

1  p-1  1  p 

Proof.  Let  Y=  (yT,...,y„)  with  y,,...,y  i.i.d.  and 
y.,  N  (0,1  >.  Then  Y/(trY’Y)^^^  =  U  =  u  AV  and  A(Y'Y/trY'Y)  = 
diag(A^, . . . ,Ap).  Note  that  A(Y’Y/trY'p  =  A(Y'Y)/tr(Y’Y)  and 
Y'Y  E  W  'P  W  (n,  I  )  ,  tr(Y’Y)  =w.,+*-*4w  and  A(W)  =  diag(w,  ,  . . .  ,w  )  , 
the  first  part  of  the  theorem  follows.  To  check  (3.5),  we  note  that 
(w^,...,Wp)  has  the  following  pdf  (cf.  Anderson  (1958)  or  Zhang  and 
Fang  (1982)): 
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(3.6) 


TT 


p/2 


2np/2  P;; 

a=0 


r(^)r(^) 


p 

n 

i=l 


l/2(n-p-l) 

w . 
i 


n 

i<  j 


(w .  -w ,  )  e 
1  3 


"2 


•+w  ) 
P 


w,  >***>w  >0. 

1  P 


Taking  the  transformation 


=  w^/ (w^+» • ‘Wp) ,  i  =  , 


X  =  Wt  +  *  •  ‘H-w  , 
1  P 


the  Jacobian  is  XP  Let  X  =  1-X.  X  Now  (3.6)  becomes 

P  1  p-1 


ttP^^  r(np/2)  P  ^%(n-p-l) 

P-X  ,._T  1  lj<i<j<p  ^  3 


n  x^ 

n  r(^)r(^) 

a=0 


n  (x.-x.) 


2np/2j,(^) 


Xnp/2-1  -X/2 


which  completes  the  proof.  Q.E.D. 


Assume  X  =  U^AV  £  (cf.  Lemma  9),  then  X  has  a  pdf 

f^(X(X'X))  if  A  has  a  pdf  f^(X^, . .  .  jX^) ,  and  f^^  is  related  to 

f^  as  follows  (cf.  Theorem  13.3.1.  of  Anderson  (1953)) 

X 


f^(X^,. . .,Xp) 


TT 


p(n+l)/2 


n  r(^)r(^) 

a=0 


(Xi  ••• 


%(n-p-l) 


n  (X . -X . )  f  (diag (X, ,...,X)),  Xt>...>X  >0. 

.  .  1  1  X  ^  ’  p  ’  1  p 

i<3  ~ 


Theorem  8.  Assume  X  e  F  with  independent  columns  (or  rows),  then 
-  ~  s 


X  must  be  normal. 
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Proof.  From  the  assumption,  the  c.f.  of  X  is  11^  d)(tlt  ),  i.e. 

~  1  ~i~i  ’  ’ 

X  E  ^2’  By  Theorem  6,  X  e  The  assertion  follows  from  Kelker 

(1970).  Q.E.D. 

This  theorem  shows  that  we,  in  general,  should  consider  dependent 

sample  theory  in  F^.  If  X  £  F^,  the  c.f.  of  X  has  the  form 

4)(X(T'T))  =  $(diag(X^,  . .  .  ,Ap))  =  tf)  (X^+*  •  •+Xp)  ,  i.e.,  the  c.f.  is  the  function 

^]_  +•  •  •+  X^.  Here  X^,  . . .  ,X^  are  the  eigenvalues  of  T'T  in  Ee^^  J  5. 

We  may  consider  other  functions  of  X, ,...,X  to  obtain  other  different 

1  P 

subclasses  of  F  . 

s  - 

5.  Applications. 

Let  F^  =  {X  e  F^:  P(X=0)=0},  ^2  ~ 

and  =  {xeF^:  P(  |x'X  |=0)=0} .  We  call  a  statistic  t(X)  distribution 
free  on  F^  if  t(X)  =  for  any  X,Y  £  F^;  i  =  1,2,3,  respectively. 


Theorem  9.  Suppose  t(X) 

is 

a  statistic.  Then 

(a) 

t(X) 

is  distribution  free 

on 

iff 

t(ap  = 

t(X) 

for 

each  a  >  0; 

(b) 

t(X) 

is  distribution  free 

on 

iff 

t(Xr)  = 

t(X) 

for 

each 

r  =  diagCr  ...  ,r  ),  r.  >0,  l  =  l,...,p;  and 

H  X 

(c)  t(X)  is  distribution  free  on  F^  iff  t(XA)  =  t(X)  for  each  A, 
an  upper  triangular  matrix  with  positive  diagonal  elements. 

Proof.  We  only  prove  (b);  the  other  proofs  are  similar.  The  "only 
if"  part  is  trivial.  Suppose  X  £  F2  and  X  =  U2R,  where 
R  =  diag (R^, . . . ,R^) ,  R^  >  0,  i  =  l,...,p,  and  U2  is  independent  of  R. 
Then  for  each  Borel  function  h  ^  0,  we  have  (by  using  the  assumption 
that  t(Xr)  =  t(X)  for  each  r) 
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E(h(t(X)))  =  E(h(t(U2R)))  =  Ej^(E(h(t(U2R))  1r)) 

=  E^(E^^(h(t(U2R))))  =  Ej^(Ey^(h(t(U2)))) 

=  E(h(t(U2)))  , 

which  is  independent  of  X  in  F^,  the  sufficiency  follows.  Q.E.D. 

Remark  1.  The  assertions  (a)  and  (c)  are  essentially  from  Kariya 
(1981a),  but  the  statement  here  is  a  little  different  from  his. 

Remark  2.  In  this  paper  we  have  only  studied  F^,  f and  F^ 

for  the  case  of  central  standard  spherical  matrices.  If  we  consider 

the  following  transformations:  F^  {x+M:  X£  F^} ,  {xE^^Vm:X£  F^} 

X  /  2 

and  F^  ^  {x^  +M:  X£  F^} ,  where  M  is  an  nXp  constant  matrix  and 

1/2  1/2 

2  =  Z  E  is  a  positive  definite  matrix,  we  can  generalize  our 
results. 

In  the  rest  of  this  section  denote  X  £  (i  =  l,2  or  3)  and 

W  =  X'D  X,  where  D  =  I  -1  1  /n  and  1  =  (1,...,!)'.  For  convenience 

~ii~  ’  ~n  ~n  ~n~n  ~n  ’  ’ 

of  applications,  some  basic  invariant  statistics  in  F^,  in  and  in 

F^  are  listed  in  Table  1.  They  are 

(1)  The  Wilks  statistic.  Let  X  =  (X' ,Xi)’,  X,:  n^  x  p,  X„:  n„xp, 

n^  >  p,  and  n^  >  p.  Let  W,  =  X/D  X,  ,  k  =  1,2.  The  Wilks  statistic 
1  —  2  —  ^  ~k  ~k-ni~-k’ 

k 

is 

t3_(p  =  |w^l/lWj^+W2l  . 

(2)  The  multivariate  Beta  statistic.  Let  W^  =  Wj^+W2,  where  W^ 

and  W„  are  given  in  (1) .  The  multivariate  Beta  statistic  is 

~  z 
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(3)  The  Hotelling  T  -statistic.  It  is 


t^CX)  =  n(n-l)x'W  , 


where  x  =  X*1  /n. 

Table  1. 

Distribution-free  Properties  of  the  Invariant  Statistics 


Statistics 

^2 

q(X) 

free 

free 

free 

tjCX) 

free 

free 

free 

t3(M 

free 

free 

free 

t4(X) 

free 

free 

free 

t5(X) 

not 

free 

free 

t^(X) 

not 

free 

free 

ft 

><1 

not 

not 

free 

rt 

00 

X 

not 

not 

not 

ft 

X 

not 

not 

not 

(4)  The  statistic  testing  equality  of  covariance  matrices. 

Partition  X=  (Xj^, .  . .  ,X^)  '  ,  where  X^:  n^^xp,  n^  p,  k=l,...,r. 

Let  W,  =  X'D  X^  ,  k=l,...,r  and  =  EW,  .  The  statistic  is 
~k  -k-n,  ~k  ~0  ~k 

k 


19 


r  n  /2  , 

t^(x)  =  n  |W^|  ^  /IWoT^^ 


(5)  The  correlation  coefficients.  It  is  easy  to  see  that  the 
sample  correlation  coefficient  between  x.  and  x.  can  be  expressed 
as 


r.,  =  xlD  x./(xlD  x.x!D  x.) 
13  ^i~n~3  ~i~n~i~j~n~3 


1/2 


Let 


t^(X)  =  R  =  (r_) 


(6)  The  canonical  correlation  coefficients.  Partition  W  into 

W  =  |„  „  I  ,  qxq,  ^22^  (p-q)  x  (p-q)  • 

The  canonical  correlation  coefficients  tg(X)  are  the  eigenvalues  of 
^12^22^21^\l‘  we  get  the  multiple  correlation  coefficient. 

(7)  Testing  the  hypothesis  that  a  covariance  matrix  is  proportional 


^11 

W 

~12 

^21 

W 

”22 

to  a  given  matrix.  Assume  Cov(x^^)  =  E  exists;  the  statistic  testing 
H:  I  =  >0  is 


t^CX)  =  lE“^wl/(tr(EQ^W)/p)P  . 


(8)  Testing  the  hypothesis  that  a  covariance  matrix  is  equal  to 
a  given  matrix.  The  statistic  testing  H:  E  =  E^  >  0  is 


tg(X)  =  [wEQ^|’^^^exp(-itr(W  Eq^))  . 


20 


(9)  The  generalized  variance. 


t  (X)  =  |X'D  X|  . 
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